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Abstract 

Control of cell proliferation is a fundamental aspect of tissue physiology central to morphogenesis, 
wound healing and cancer. Although many of the molecular genetic factors arc now known, the system level 
regulation of growth is still poorly understood. A simple form of inhibition of cell proliferation is encountered 
in vitro in normally differentiating epithelial cell cultures and is known as " contact inhibition" . The study 
presented here provides a quantitative characterization of contact inhibition dynamics on tissue-wide and 
single cell levels. Using long-term tracking of cultured MDCK cells we demonstrate that inhibition of cell 
division in a confluent monolayer follows inhibition of cell motility and sets in when mechanical constraint on 
local expansion causes divisions to reduce cell area. We quantify cell motility and cell cycle statistics in the 
low density confluent regime and their change across the transition to epithelial morphology which occurs 
with increasing cell density. We then study the dynamics of cell area distribution arising through reductive 
division, determine the average mitotic rate as a function of cell size and demonstrate that complete arrest 
of mitosis occurs when cell area falls below a critical value. We also present a simple computational 
model of growth mechanics which captures all aspects of the observed behavior. Our measurements and 
analysis show that contact inhibition is a consequence of mechanical interaction and constraint rather than 
intcrfacial contact alone, and define quantitative phenotypes that can guide future studies of molecular 
mechanisms underlying contact inhibition. 
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1 Introduction 

The precise orchestration of cell division and growth 

is central to morphogenesis and animal develop- 
ment [1, 2]. Complex cellular signaling and regu- 
latory networks are dedicated to growth control and 
misregulation of cell proliferation leads to tumors and 
cancer [3] . Epithelial tissue is an important system to 
study regulation of growth. Normal development of 
epithelial tissue involves a mesenchymal to epithelial 
transition (MET) [4] associated with the loss of cell 
mobility, mitotic arrest and acquisition of epithelial 
morphology. This transition is reversed in the process 
of wound healing [5]. On the other hand, cells that 
have undergone oncogenic epithelial to mesenchymal 
transition (EMT) typically lose their ability to un- 
dergo MET. Hence understanding the normal MET 
process is of fundamental importance for understand- 
ing oncogenic transformations which disregulate it. 

In cultured, non-cancerous epithelial cells, the 
transition from freely proliferating, non-confluent 
cells to fully differentiated, dense epithelial mono- 
layers is commonly referred to as "contact inhibi- 
tion" [6, 7, 8, 9]. Contact inhibition in confluent 
cell cultures is currently defined as i) a dramatic 
decrease of cell mobility and mitotic rate with in- 
creasing cell density; ii) establishment of a stationary 
post-confluent state which is insensitive to nutrient 
renewal. It is widely believed that contact inhibition, 
as the name suggests, is caused by cell contact. But 
despite extensive study, current understanding of the 
mechanism of contact inhibition is far from complete 
(see [10, 11, 12, 13, 14]). 

Many molecular mechanisms have been proposed 
to contribute to contact inhibition. It is widely ac- 
cepted that contact inhibition requires establishment 
of E-cadherin mediated cell-cell contacts and subse- 
quent maturation of the adherens junctions (AJs) 
that link E-cadherin and F-actin in a synapse-like 
complex involving numerous other proteins [15, 16, 
17, 18]. However, the nature of the signaling path- 
way leading to suppression of mitosis remains un- 
clear. One possible pathway involves /3-catenin, 
a mediator of Wnt signaling, that, in addition to 
its function as a transcriptional co-factor, is asso- 



ciated with the AJs at the cell surface [19, 20]. 
A contact inhibition role has been reported for 
NF2/Merlin, a tumor suppressor gene [21, 22] that 
encodes a mcmbranc-cytoskeletal scaffolding protein, 
which most likely acts via the Hippo kinase path- 
way, controlling nuclear localization of the transcrip- 
tional activator YAP [23, 24, 13] - itself a known 
regulator of cell proliferation. Contact inhibition 
is known to involve the MAPK pathway, which, in 
turn, promotes cell cycle entry by regulating the ex- 
pression of cyclinDl [25, 26, 27]. Also implicated 
are Nectins [28, 29, 12] - a family of cell adhesion 
molecules that are involved, together with integrins 
and other proteins, in the regulation of cell motility 
and proliferation. Yet, this accumulated knowledge 
falls far short of a comprehensive picture of contact 
inhibition. The difflculty in achieving a better under- 
standing of the molecular mechanism lies in the com- 
plexity of the contact inhibition phcnotypc, which, as 
we describe below, involves the concurrence of many 
processes. 

To facilitate progress in the dissection of the reg- 
ulatory pathways underlying contact inhibition, we 
undertook a quantitative reexamination of the spatio- 
temporal dynamics of an adherent epithelial layer 
formed by Madin-Darby Canine Kidney (MDCK) 
cells. These cells are known for their ability to exhibit 
contact inhibition and achieve characteristic epithe- 
lial morphology in culture [30, 31] thus providing an 
excellent model system for in vitro study of epithelial 
tissue dynamics [32, 33, 34]. Using long-term fluo- 
rescence and phase contrast video microscopy in con- 
junction with image segmentation and cell tracking, 
we have characterized the temporal progression of 
contact inhibition in growing MDCK colonies. Quan- 
titative analysis of the evolution of cell density, cell 
motility and cell division rate reveals that contact 
inhibition proceeds in three distinct stages: 1) a 
stage of cell density growth with gradual inhibition of 
motility, but without inhibition of mitosis that is fol- 
lowed by 2) a rapid transition to epithelial cell mor- 
phology, followed by 3) continued cell division and 
reduction of cell size with a progressively decreasing 
rate of mitosis. Mitotic arrest is achieved once cell 
area falls below a certain threshold. Our findings 
show that contact between cells is not sufficient for 
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inhibition of mitosis in MDCK cells. Instead, inhibi- 
tion of cell proliferation is a consequence of mechan- 
ical constraint that causes successive cell divisions to 
reduce cell area. 

2 Results 

2.1 Large scale analysis 

To separate the effect of cell contact from that of me- 
chanical constraint arising upon confluence of prolif- 
erating cells, we first examine the dynamics of iso- 
lated, growing colonies of MDCK cells. The colonies, 
started from a small initial number of cells, were mon- 
itored with subcellular resolution by time-lapse video 
microscopy for up to three weeks until nearly com- 
plete proliferation arrest (see Materials and Methods 
for details). Fig. lA shows the large scale dynam- 
ics of a growing colony. The boundary of the colony 
exhibits nontrivial dynamics due to the combined ef- 
fect of motility and cell division. It moves outward 
with a non-uniform velocity forming finger-like pro- 
trusions [32, 35, 33]. Yet the total area of the colony 
grows following a simple exponential law (Fig. IB) 
for up to 5-6 days, reaching over 10'^ cells. Cell den- 
sity in the bulk remains constant during this period 
(Fig. IC). Daughter cells occupy, on average, twice 
the area of their mother cell and the rate of colony 
area growth matches exactly the rate of cell mito- 
sis. Thus, colony expansion is driven by cell prolif- 
eration. In this "free growth" regime, although cells 
stay in contact with each other for several days, no 
inhibition of growth is observed. 

Expansion of the colony is made possible by the 
fluid-like motion of cells in the two-dimensional con- 
fluent layer (see Fig.2A). To quantitatively charac- 
terize this motion we carried out a PIV-type analysis 
of the time series of phase-contrast images [36] . This 
analysis determines the local velocity field by compar- 
ing successive images. In the free growth regime, cells 
exhibit a swirling, but outwardly biased flow with a 
root-mean-square (r.m.s.) outward velocity of about 
15/xm/h, or approximately 1 cell width per hour. The 
motion of nearby cells is correlated on a length scale 
of about 5 cells (see Fig.2C). 



Of course, exponential increase in colony area can- 
not continue indefinitely as it would require the out- 
ward motion of peripheral cells to have an exponen- 
tially increasing velocity. To support an exponential 
area increase at a rate 1/r, the velocity of cells on 
the boundary must be Vb — \fAj^KT^ . Comparing 
this to the observed maximal velocity of cell motion, 
■f^max — 15/zm/h, we arrive at an estimate for the crit- 
ical size of the colony, Ac 2 ■ Above this 
critical area, expansion of the colony cannot keep up 
with cell proliferation in the bulk without increasing 
cell density. This estimate is close to the observed 
area at the time of crossover (~ 6 days) from expo- 
nential to sub-exponential growth of the colony area 
(Fig. IB). The crossover is indeed coincident with the 
onset of a gradual increase in cell density in the bulk 
of the colony, as shown in Fig. IC. Single cell anal- 
ysis (below) confirms that mitotic rate in this "pre- 
transition" regime does not decrease so that the sub- 
exponential expansion of the colony area is accounted 
for by increasing density alone. 

Thus, it appears that increased cell density (and 
the associated decrease of average cell size) is a con- 
sequence of mechanical constraint imposed by the in- 
ability of the tissiie at the periphery to expand fast 
enough to accommodate cell proliferation in the bulk. 
As cell density begins to increase, cell motility starts 
to decrease as shown in Fig. 2C (see also ref. [33]). 
The correlation length associated with the velocity 
field exhibits a peak which closely corresponds to 
the transition to epithelial cell morphology, which 
we shall define and discuss in detail below. The 
correlation length of motion subsequently decreases 
with time down to the size of a single cell (where the 
displacement is comparable to the optical resolution 
of the images), indicating that large scale swirls ob- 
served in the free growth regime are disappearing and 
cell motion is reducing to small scale fluctuations. 
The observed decrease in the root-mean-square ve- 
locity and the correlation length of cell motion can 
be understood in terms of a rapid increase in cell- 
substrate adhesion in the process of morphological 
transition. 

Cell behavior in the center of a colony at the end of 
the free growth regime is similar to what is observed 
in confluent cell cultures that were seeded homoge- 
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ncously. sec Fig.2C,D.In homogeneously seeded cul- 
tures the space constraint is more severe and cell den- 
sity increases more rapidly upon confluence, reaching 
the morphological transition soon after initiation of 
the culture. We note that panels C and D in Fig. 2 
differ also in the substrate: glass and PDMS respec- 
tively (see Materials and Methods). Although the 
magnitude of the correlation length differs in the two 
cases, the qualitative behavior and the morphological 
transition are very similar. 

To investigate the long term evolution of cell size 
and density upon exit from the free growth regime we 
measured the average cell area (using digital image 
segmentation) over a period of 20 days. As shown in 
Fig.SC, average cell area decreases 10-fold over a 15 
day period. In the same period, average cell height 
increases only by a factor of 2 (from 5-6 /xm to 12-15 
/Ltm, see Supplementary Information Fig.SlD) thus 
indicating that cell volume decreases. Fig. 3C shows 
a rather sharp crossover from rapid to slow growth of 
cell density (and corresponding decrease of cell size). 
This crossover coincides with the transition in cell 
morphology illustrated by the comparison of panels 
A and B in Fig. 3 and separates the "pre-transition" 
and "post-transition" stages of the contact inhibition 
process. The pre-transition transient is characterized 
by the gradual reduction of cell motility (Fig. 2) 
discussed above. The post-transition state is charac- 
terized by the absence of cell rearrangement, except 
through cell division. Mitotic rate decreases contin- 
uously in the post-transition stage (see Single cell 
analysis section) leading to arrest of cell prolifera- 
tion independent of nutrient renewal. This "arrested" 
regime can last for weeks. (We monitored the tissue 
for 23 days without detecting any significant changes 
in the area distribution.) However, the state of pro- 
liferation arrest can be readily reversed by scratching 
the cell layer to create a free boundary [32], or by 
stretching the substrate along with the cell layer. 

The morphological transition itself is readily quan- 
tified by the radial distribution function which mea- 
sures conditional probability of finding a cell at a 
given distance from a reference cell (Fig. 3D). In 
the pre-transition stage, the distribution function ex- 
hibits an exclusion zone at distances comparable to 
the size of the nucleus, and is flat for larger distances. 



indicating a disordered system of non-overlapping 
cells. In the post-transition regime a peak and a 
trough emerge in the distribution function, corre- 
sponding to nearest and next-nearest neighbors. This 
indicates an increase in size homogeneity and the 
appearance of local ordering of cells within the tis- 
sue. Further quantitative characterization of the 
these regimes is provided by single cell tracking and 
analysis. 

2.2 Single cell analysis 

To further quantify cell behavior in the different 
regimes of tissue growth, we followed individual cells 
through the division process. In the free growth 
regime, each daughter cell grows back to the size of its 
mother cell, as shown in Fig. 4A. In contrast, in both 
pre- and post- transition regimes mitosis reduces cell 
area by approximately a factor of two without sub- 
sequent growth of daughter cells (Fig 4A, B). The 
prevalence of reductive division is demonstrated in 
Fig. 4C, which compares measured areas of mother 
and daughter cells. Combined area of daughter cells 
does not exceed the area of the "mother" cell, in- 
dependent of the area of the latter (Fig. 4C, in- 
set). Thus, cells in both the pre- and post-transition 
regimes are "incompressible" in the sense that new 
cells introduced through cell division do not achieve 
any expansion of their area at the expense of their 
neighbors. 

The morphological transition coincides with an ap- 
proximately 5-fold decrease in the average mitotic 
rate and a dramatic broadening of the distribution 
of cell cycle periods (Fig. 4E). Once the cell divi- 
sion time becomes significantly longer than a day, 
measuring mitotic rate by tracking individual cells 
becomes very difficult, necessitating a different strat- 
egy for measuring the dependence of mitotic rate on 
cell size. The latter, as we show next, can be deduced 
from a quantitative study of the temporal evolution 
of the cell area distribution. 

To measure the distribution of cell area as a func- 
tion of time in the post-transition regime, we per- 
formed a computerized segmentation of fluorescent 
images (using the MDCK-EcadGFP ceU line) (Fig. 
3B and Fig. 4D). Over a period of about 15 days fol- 
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lowing the morphological transition the average cell 
area decreases 6-fold. Cell size converges on a narrow, 
stationary size distribution centered about an aver- 
age area of 35^m^ (see SI, Fig. SI for more details on 
cell morphology). 

Because cell area in the post-transition phase is 
approximately constant between successive divisions 
(Fig. 4A-B), the dynamics of the cell area distribu- 
tion is due solely to mitosis. (Rates of apoptosis, 
about 0.02 per day per cell, arc negligible by com- 
parison with mitosis.) Thus the difference in area 
distributions at two consecutive time points reflects 
loss of larger cells which upon division give rise each 
to a pair of cells at approximately half the size as 
represented by: 

dtn{a, t) = 2j{2a)n{2a, t) - 7(a)n(a, t) (1) 

where n{a,t) denotes the expected number of cells 
with area, a, at time t and 7(a) represents the average 
rate of division as a function of cell area. Fitting the 
observed temporal changes in cell area distribution to 
Eq. (1) allows us to determine the mitotic rate 7(a). 
The dependence of mitotic rate on cell size is shown in 
Fig. 4F. The result is consistent with the conclusion 
made on the basis of the single cell measurements, 
Fig. 4E: a rapid decrease in the mitotic rate once cell 
area falls below critical, which fits approximately the 
Hill function form: 7(a) /70 = a^^a" + ag") with 
m « 4 and ao « 170/xm^, where 70 is the division 
rate in the free growth regime. The observed reduc- 
tive nature of cell division and the dependence of mi- 
totic rate on cell size together explain the dynamics 
of tissue density in the post-transition regime and the 
convergence towards proliferation arrest, thus captur- 
ing the dynamics of the "contact inhibition" process. 

2.3 A model of self-limiting growth of 
adherent cell monolayer 

To illustrate our interpretation of the observed in- 
terplay between cellular growth, motility and colony 
expansion, we formulate and analyze a simple growth 
model for adherent epithelial tissues. We choose as a 
point of departure a one dimensional version of the 
"vertex model" [47, 50] as depicted in Fig. SA.The 



details of the model arc described in the supplemen- 
tal material. Briefly, we assume that cells, specified 
by their vertices and rj+i, form a connected tissue. 
The short time clastic response of cells has a Hookc's 
law dependence (see SI Eq 2.1) on the difference be- 
tween the current length of the cell, k = |rj+i — r,! 
and the intrinsic preferred length L; . In addition cells 
interact with the substrate. To represent the effect of 
substrate adhesion and of cell motility we introduce 
for each cell an "attachment point" Ri, connected to 
the cell by a spring, and endowed with relaxational 
dynamics with friction a and random Langevin driv- 
ing force rii{t) (see SI Eq 2.2). This (gaussian white) 
random force represents cell crawling, its variance T 
defining motility. The force is assumed to average out 
to zero in the bulk, but not on the boundary where 
to represent the outward bias of the boundary cell 
motion we allow < rjb >= crvmax, where Vmax sets 
the maximal velocity. 

The model also includes cell growth and prolifera- 
tion. Cell growth is represented by allowing intrinsic 
cell size Li to increase with time. However, moti- 
vated by our experimental finding that cells in dense 
tissues don't actively push on their neighbors, we al- 
low Li to grow only if the cell is stretched by the sur- 
rounding tissue. Because stretching corresponds to 
Sli = li — Li > -we take dLi/dt to be a simple step 
function of Sk with the threshold at zero. Cell divi- 
sion splits a cell into two, with intrinsic size of each 
daughter equal to Lj/2 of the mother. Guided by our 
experimental observations (Fig. 4) we make the rate 
of cell division explicitly dependent on cell size li via 
p{li) = max(7(?j — Zinin);0) which implements a "size 
check point" (li > /min) for cell proliferation. 

The spatio-temporal growth dynamics of an ini- 
tially small colony is shown in Fig. 5B. Independent 
of the specific parameters used in our simulations, 
we find two growth regimes: an initial exponential 
growth of the colony with uniform proliferation is 
followed by a linear growth regime with cell prolif- 
eration mainly on the margin of the colony. The 
maximal size of a freely proliferating colony and the 
critical time for cross over from exponential to sub- 
exponential growth may be obtained by the following 
argument. A patch consisting of N{t) = e^* prolif- 
erating cells increases its diameter with the instan- 
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tancous speed lodN{t)/dt. This speed cannot exceed 
the maximal velocity of the interface Vmaxj lest pro- 
liferating cells become compressed, as the colony be- 
comes unable to spread fast enough to keep up with 
cell proliferation. The two velocities became equal at 
tc = iog{2vmax/'yio) 8kt which time the colony size 
is 2vmaxh- (This ID argument is readily generalized 
to 2D.) 

The two growth regimes are manifest on the colony 
scale but are also reflected on a single cell level (Fig. 
5B inset). Large cells on colony margins re-grow in 
size after division and then divide again. In higher 
density regions, cells undergo size reductive divisions 
until the division finally ceases. The temporal evo- 
lution of the cell sizes in the middle of the colony 
is shown in Fig. 5C. The diS'crence in cell behavior 
arises, in our model, from the difference in their me- 
chanical state. Motile cells in small colonies or on 
the periphery of a large colony arc under tension and 
grow in size after division. In contrast, cells in the 
bulk of a large colony are "boxed in" by their neigh- 
bors, and after each division reduce in size until they 
divide no more. Our assumptions that a) cells grow 
in size only under tension and b) cells do not exert 
compressive forces on each other, are closely related 
and result in a distribution of tensile stress across the 
tissue layer (see Figure 5D) that is consistent with the 
observations of Trepat et al [45]. 

3 Discussion 

Quantitative observations of cell size, motion and di- 
vision rate reported here help dissect the complex 
nature of the "contact inhibition" phenotype. They 
reveal, for example, that free exponential growth can 
take place within cell colonies even after cells have 
been in contact for several days. Such a long delay 
between confluence and mitotic inhibition is partic- 
ular to isolated, expanding colonies. Onset of mi- 
totic inhibition occurs sooner in cell cultures seeded 
at uniform density, where confluence coincides with 
near complete occupation of the available area. We 
conclude that cell-cell contact is a necessary [37] but 
not sufficient condition for growth inhibition. 

The data strongly suggests that inhibition of cell 



division follows the reduction in cell area imposed by 
mechanical constraints on tissue expansion [38, 37]. 
Interestingly, although average cell area starts to de- 
crease in the motile prc-transition regime, the depen- 
dence of mitotic rate on cell area does not appear 
until the cells enter the static post-transition phase, 
at which point the rate of cell proliferation drops 
sharply. The cessation of cell motility and drop in 
mitotic rate coincide with a dramatic change in tis- 
sue morphology, as revealed by the radial distribution 
function of cells. The radial distribution function is 
thus an interesting, quantitative characterization of 
the state of the tissue and its abrupt change can serve 
to pin-point the morphological transition, which may 
be considered a key element of the mesenchymal-to- 
epithelial transition (MET). 

Our study focused on the dependence of mitotic 
rate on cell area because the latter is directly measur- 
able. Yet it is important to emphasize that the real 
trigger of intracellular signals responsible for the sup- 
pression of motility, inhibition of mitosis, and MET 
may be not cell size per se, but mechanical stress 
and deformation which are known to induce reorga- 
nization of focal adhesion, adherens junctions, and 
cytoskeleton [10, 39, 40, 41, 42]. 

It is widely believed that mechanical tension pro- 
motes cell division [43, 44]. Recent direct measure- 
ments by Trepat and coworkers [45] demonstrate that 
motile MDCK cells exert inward pointing traction on 
the substrate, implying that spreading cell colonies 
are under tension. These measurements correspond 
to MDCK colonies in what we refer to as the free 
growth regime. Combining their observations with 
ours suggests that the motile pre-transition regime 
and the morphological transition itself correspond to 
the gradual relief of tension and the onset of com- 
pression brought on by cell proliferation. It has 
been suggested (in the context of the problem of 
organ size determination) that mechanical compres- 
sion may be providing an inhibitory signal for mito- 
sis [46, 47, 48, 49]. The same hypothesis could then 
explain the observed inhibition of mitosis in the post- 
transition regime. If so, the trigger of MET would be 
the change from tensile to compressive stress acting 
locally within the cell layer. Our ID model illustrates 
how this scenario can generate the observed behavior. 
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Our measurements also suggest that inhibition of 
cell division is a distinct single cell state rather than 
a global state induced by cell-cell signaling across 
the layer, as illustrated in Fig. 4E. In fact, confluent 
MDCK cell cultures with an average cell density cor- 
responding to the morphological transition are often 
sufficiently heterogeneous in local cell density that 
highly motile cells and completely arrested cells co- 
exist in the same colony. Thus contact inhibition is a 
local phenomenon, which calls into question the reli- 
ability of "bulk assays" of the phenotype. Hence fu- 
ture experimental efforts focusing on molecular mech- 
anisms underlying mitotic inhibition should be con- 
ducted with methods allowing single cell resolution. 
Such experiments, combined with techniques allow- 
ing in situ measurements of mechanical stress acting 
on cells would, we believe, finally lead to decisive un- 
derstanding of the contact inhibition phenomenon. 



were identified by the fully automatic segmentation 
of Ecad-GFP fluorescence images. Mean displace- 
ment measurements were made by cell tracking using 
a PIV-type analysis [36]. Colony profiles were ob- 
tained by extracting the boundary of the colony as a 
function of time from the images using standard edge 
detection algorithms. 

Numerical simulations The mathematical model 
underlying Fig. 5 and the method of simulation are 
described in the SI. 
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4 Materials and Methods 

Cell culture MDCK-II cells and MDCK-Ecad-GFP 
were a gift from J. Nelson. Cells were cultured 

in MEM (GIBCO, 11095-098) supplemented with 
Penicillin-Streptomycin and 5% FES (Cellgro, 35- 
010-CV) at 37° C and 5% CO2. 

Time-lapse microscopy All imaging was per- 
formed on an inverted microscope (Olympus IX- 70) 
with a 20X/0.7NA/Ph2 objective. Phase and fluores- 
cence images were taken with, respectively, a halo- 
gen lamp and an LED (Luxeon LXHL-LB5C) and 
recorded to disk using a CCD camera (Qlmaging, 
Retiga EXi). Mechanical shutters (Uniblitz VS25) in 
both illumination paths limited sample exposure to 
minimize phototoxicity. Stepper motors controlled 
the stage position and objective focus. The shut- 
ters, stepper motors and camera were controlled by a 
custom-written Labview program. Images were taken 
every 10 min (phase) and every 3h (fluorescence) for 
a given field of view. 

Image analysis Images were analyzed by using cus- 
tom written Matlab programs. Positions of nu- 
clei in the low cell density phase were determined 
from the phase contrast images, segmented semi- 
automatically. In the static, high density regime, cells 
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Figure 1: Epithelial colony growth. A) Superimposed snapshots of a single colony at different times, 
coded by different shades of gray. Time-points were chosen to keep area increment constant. Black contours 
correspond to 3.0, 4.8, 5.5, 5.9, 6.3 days after seeding. B) Total area of the spreading colony. Time is 
counted relative to the "morphological transition" at i = (see text). Green points represent total cell 
number (independently measured) multiplied by the average cell area. The blue line is exponential growth 
with the average cell cycle time T2 = 0.75 ± 0.14 (s.e.m.) days (measured independently by single cell 
tracking). C) Cell density in the inner region of the colony (different colors distinguish different fields of 
view). The solid black line at constant density and is a guide for the eye. The dashed black line represents 
exponential growth of density expected for continued cell proliferation without cell motion. 
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Figure 2: Correlation analysis of cell motility. A-B) Phase images of a confluent layer (Ih before 
and 27 h after the morphological transition) with overlaid instantaneous velocity field (measured by PIV 
and interpolated) side by side with cell trajectories integrated over 200 minutes with blue and red labeling 
respectively the beginning and the end. Scale bar is lOO/im. C-D) R.m.s. velocity of cell motion (red 
symbols) and the correlation length (blue symbols) across the morphological transition in a in the bulk of 
the expanding colony (C) and in the continuous confluent layer plated at higher initial density (D). Data 
pooled from four different 450x336 fj,m? fields of view. Lines are to guide the eye. D) inset: Correlation 
time of cell trajectories. 
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Figure 3: Large scale quantitative characterization of contact inhibition. A- B) Image segmentation 
for MDCK cell cultures grown on PDMS. A) phase contrast image of Ecad-GFP MDCK at low cell density. 
B) fluorescent image of Ecad-GFP MDCK at high cell density. Scale bar: 100/im. C) Median of cell 
area distribution (over a 450x336 iirri^ field of view) as a function of time (t=0 set to the morphological 
transition). Here MDCK cells were seeded at uniform density (see methods) and imaging commenced 
upon confluency. Different line colors represent different experiments which are time aligned. No density 
change is detectible after 15-18 days. D) Radial distribution function of cells at different times across the 
morphological transition. g{R) is the ratio between the density of cells in a circular annulus distance R 
from a reference cell and the average density. The appearance of a peak (and a trough) in the static post- 
transitional phase represents increased short range ordering of cells. t=0 is defined by the first appearance 
of the peak: max{g{R)) > 1.2. 
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Figure 4: Single cell level quantification of contact inhibition. A) and B) Traces of single cell area 
tracked as a function of time. Arrows represent mitosis. A) starts below confluence and reaches high density 
confluence; B) is in the post-transition phase. Dashed lines represent temporal averages. C) Daughter cell 
area versus the area of the mother cell. Data represent 96 divisions at different times for confluent layers. 
. Daughter cell area as the average over three time points 1 hour apart, 12 hours after mitosis. Mother vs 
daughter cell areas follow the line y = x/2, plotted in black. C) Upper inset: Distribution of daughter 
cell areas in % of the mother cell area. C) Lower inset: Deviation of the total daughter cell area from 
the mother cell area. D) Distribution of cell area in the post-transition regime. Color codes for time. Each 
distribution represents the population of (at least 200) cells in the same 336x450 /im^ field of view. Cell area 
is measured by means of computer segmentation of MDCK-Ecad-GFP fluorescent images (see fig. 5A). E) 
Single cell division times as a function of premitotic area. Different colors represent different experiments. 
We note that cell cycle time increases dramatically for cell areas below 200 /xm^ . The absence of data below 
lOfim^ is due to the difficulty of tracking single cells in that regime. E) inset: Distribution of division times 
in the pre-transition regime. F) Division times as a function of cell area inferred from the dynamics of the 
P{A) functions using Eq. (1) (see methods). Different colors represent different experiments. The black line 
represents average division time in the pre-transition regime. Cell division slows down for cell size below 
200/Ltm^, consistent with cell tracking measurements shown in panel E. 
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Figure 5: Simulation results. A) Sketch of the one-dimensional tissue growth model. Green springs 
represent cell elasticity, cell boundaries are marked in gray and cell attachments are represented as black 
tethers. B) Spatio-temporal profile of proliferation rate (indicated by the color) in the colony. Initially, 
proliferation is uniform and the colony size increases exponentially (dashed line) with time. At later times, 
proliferation in the bulk slows down and stops; in fixed size marginal zones rapid cell proliferation continues 
leading to a linear increase of the colony size (dotted line). C) Cell size distribution as a function of 
time (coded by color). The initial distribution around (set by the ratio of the rates of cell growth and 
proliferation) becomes broader and converges with time to a stationary distribution with mean below Imin- 
Inset shows the coefficient of variation. D) Traction force distribution throughout the colony at different 
times (coded by color) . Note that small colonies are under tension. At later times only the margins of the 
colony are under tension, while the center is stress free. 
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1 Supplementary data 



Movie SI.l: Dynamics of a growing epithelial colony. The movie is composed of 
9x9 patches for a total of 3,02x4,05 mm. 6 cells were seeded and imaged for 10 
days. Media were replaced daily. Scale bar: 1mm. 



Movie SI.2: Dynamics of epithelial tissue in the bulk. Cells were seeded at 
uniform density and imaged continuosly for a week. Media were replaced daily. 
Confluency is reached at t = 0.6 days. By day 2 cell motion has completely 
disappeared, and only small scale vibrations can be observed. Scale bar: lOO/xm. 
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Figure SI.l: A) Aspect ratio of the Voronoi cells constructed for the segmented 
nuclei. Aspect ratio is defined as the ratio between the largest and the smallest 
eigenvalues of the inertia tensor defined by cell vertices. Regular polygons have 
aspect ratio between f and 1.5. The figure presents the fraction of the cell 
population with aspect ratio ar < 1.5 (black), 1.5 < ar < 3 (blue) and ar > 3 
(red). As a function of time, the fraction of regular polygons increases and 
the fraction of deformed cells decreases dramatically. Time zero refers to the 
morphological transition. B) Change in cellular the coordination number across 
the morphological transition. Each histogram (coded by a color) represents the 
distribution of the number of cell vertices at a given time during the experiment. 
The last distribution plotted is stationary and is analogous to the one found 
in other biological tissues [3, 4, 5, 6]. C) Cell size heterogeneity across the 
morphological transition. Different colors represent the coefficient of variation 
as a function of time for different experiments. The coefficient of variation was 
obtained by calculating the ratio of the difference of the quartiles 75 and 25 
to the median: (Q75 — Q2b)IQ50- Data from different experiments were time 
aligned referring to the morphological transition. D) Average cell thickness as 
a function of cell area. We observe that below 200/im, cell height does not 
change appreciably. Cell height is measured by finding the best focal plane at 
a given time and then tracking the position of the objective during the whole 
experiment. 
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2 Langevin model for tissue motility. 



To illustrate our interpretation of the observed interplay between cellular growth, 
motility and colony expansion, we formulate and analyze a simple growth model 
for adherent epithelial tissues. We choose as a point of departure a one dimen- 
sional version of the "vertex model" as shown in Fig. 5A. Physical position of 
cell i is specified by its vertices and r^+i which it shares with its neighbors. In 
addition to the interaction between neighboring cells we also include adhesion 
to the substrate. The attachment of cells to the extracellular matrix is medi- 
ated by many focal adhesions on the basal membrane of the cell. In our model 
we represent the cumulated effect of the focal adhesions by a single attachment 
point Ri for each cell. We assume the mechanical properties of the tissue to be 
elastic on short time scales so that for a given set of intrinsic cell lengths and 
attachment points, the vertices of the cells are determined by minimizing the 
energy 

N N 

H{ri, VN+i) = k ^(n+i -n- Li{t)f + K - in+i + r•i)/2]^ (SI.l) 

i=l i=l 

where Li{t) is the intrinsic preferred length of the cell. The first term describes 
the mechanical interaction between cells, the second term accounts for the at- 
tachment of the cells to the substrate and the units are chosen so as to make 
attachment stiffness unity. In principle, depending on the mechanical stress 
Si oc k{ri^i — Ti — Li{t)) a cell can adapt its intrinsic length scale Lj resulting 
in an effective plasticity of the tissue that would relax the stress at long times. 
Here, for simplicity, we will only include the effect of cell growth, which also 
manifests itself as a change - an increase - of Li. We assume that Li resists 
compression s, < and grows with rate a for Sj > 0, e.g. ^Li = a9{si), where 
6 is the step function. This implements the assumption that cell grow in size 
only when under tension, with the consequence in the absence of tension, cell 
divisions will reduce intrinsic cell size Lj. 

The dynamics of the attachment points is driven by relaxation of elastic 
stress and by random forces generating cell motility 

where rji - a Langevin-type random force representing motility - is a Gaussian, 
white random function of time defined by its second moment < r]i{t)r]j{0) >= 
r6{t)Sij. The relevant time scale for the attachment dynamics is given by cr - 
which acts like friction. While we assume that the random force representing cell 
motility has a zero average in the bulk of the tissue, to represent the observed 
directed (outward) crawling of cells on the boundary, we allow < tjn >= — < 
Tji >= crvmax, where v^ax sets the maximal crawling velocity of boundary cells. 

In addition to the continuous changes in cell sizes and attachments, cells may 
divide. Motivated by our experimental quantification of the cell area growth 
curve (Fig. 4) we make the average rate for cell division explicitly dependent on 
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cell size li = Tj-i-i — Tj by using p{li) = Taax{"/{li — ^inin),0) in our simulations, 
which is also consistent with our finding of a growth size check point. The 
division process replaces the dividing cell by two equivalent daughter cells with 
attachment points set to the middle of each cell and the sum of internal length 
matches the internal length of the mother cell. 

We have simulated the dynamics of this model (with parameters k = 1,k = 
.5. a = .02,(7 = l,r = 1,7 = .2) using custom-written Matlab programs imple- 
menting matrix inversion for the dynamics of the vertices and the Runge-Kutta 
method for the dynamics of the attachment-points. Programs are available upon 
request. 

3 Additional information on the methods 

Cell culture Cells were seeded at uniform density (around 600 cells/mm^) on a 
fibronectin (Sigma- Aldrich, F1141-2MG) coated PDMS membrane (McMaster- 
Carr, 87315K62) and imaged in phenol red free IMEM (Cellgro, 10-26-CV) sup- 
plemented with Penicillin-Streptomycin and 5% PBS. The media was replaced 
daily and the culture conditions were kept at 37° C and 5% CO2 by means of a 
custom made microscope stage enclosure. Single colony experiments were per- 
formed by seeding cells at a density of about 1 cell/cm^ in a glass bottom petri 
dish. Images of the colony spanning 9x9 contiguous fields of view were captured 
and stitched together, making it possible to follow single cell motion as well as 
tissue wide dynamics with good resolution. 

Image processing. Quantitative data on cell area was obtained with the 
help of contrast enhancement by Gabor-filtering [7, 8] and successive removal of 
low contrast regions. Frames taken 10 minutes apart were compared to remove 
poorly segmented cells. Cell size, shape, topology and structure functions were 
obtained from the segmented images. Mean displacement measurements were 
made by cell tracking using a PIV-typc analysis [9]. The position of nuclei 
in a segmented image was propagated in time by correlation analysis. The 
trajectory thus obtained was then corrected for stage drift and rms velocities 
were calculated as (< > — < w >^)^/^. To extract the correlation length, 
the two point correlation function Cv{R) for {?— < v > was computed and the 
length for which Cy{R) = 0.3 was used. 

Rate analysis To calculate the division rates shown in fig. 4E, cq. 1 was dis- 
cretized; an over-determined linear system was constructed based on measure- 
ments of the size distribution at different times, and was subsequently pseudo- 
inverted [10] . The total number of cells per frame was obtained from the average 
area. 
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